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ABSTRACT 

We report the discovery of a powerful molecular wind from the nucleus of the non-interacting nearby 
SO field galaxy NGC 1266. The single-dish CO profile exhibits emission to ±400 km s~^ and requires a 
nested Gaussian fit to be properly described. Interferometric observations reveal a massive, centrally- 
concentrated molecular component with a mass of 1.1 x 10^ Mq and a molecular outflow with a 
molecular mass of w 2.4 x 10^ M©. The molecular gas close to the systemic velocity consists of a 
rotating, compact nucleus with a mass of about 4.1 x 10^ Mq within a radius of « 60 pc. This compact 
molecular nucleus has a surface density of ~ 2.7 x lO"' Mq pc~^, more than two orders of magnitude 
larger than that of giant molecular clouds in the disk of the Milky Way, and it appears to sit on the 
Kennicutt-Schmidt relation despite its extreme kinematics and energetic activity. We interpret this 
nucleus as a disk that confines the outflowing wind. A mass outfiow rate of w 13 Mq yr~^ leads to a 
depiction timescale of < 85 Myr. The star formation in NGC 1266 is insufficient to drive the outflow, 
and thus it is likely driven by the active galactic nucleus (AGN). The concentration of the majority 
of the molecular gas in the central 100 pc requires an extraordinary loss of angular momentum, but 
no obvious companion or interacting galaxy is present to enable the transfer. NGC 1266 is the first 
known outfiowing molecular system that does not show any evidence of a recent interaction. 
Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — ISM: jets and outflows 
— galaxies: ISM — galaxies: nuclei — galaxies: elliptical and lenticular, cD — 
ISM: kinematics and dynamics — galaxies: individual (NGC 1266) 
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1. INTRODUCTION 

Early-type galaxies (ETGs) are ge nerally thought t o 
be poor in atomic and molecular gas (jLees et al.lll991[ ). 
but their disk g alaxy progenitors are generally gas-rich 
(jKauffmann et al. 2003). How these galaxies transform 
themselves from gas-rich to gas-poor systems is tied to 
a host of issues in galaxy evolution, including the migra- 
tion of galaxies fro m the blue cloud to the red sequence 
(iFaber et al.|[2007l). merger-d riven starbursts and quasar 
activity ( Hopkins et al.|[2005[) , feedback from star forma- 
tion and central supermassivc black holes on the inter- 
stellar medium fISM) (fPi Mattco et al...2005i) and radio 
mode activity (jCroton et aLll20d6l ). 
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The ATLAS^'^ project aims to address these issues 
through a complete volume-limited multi-wavelength 
survey of ETGs within 42 Mpc, spanning a variety 
of environments and tw o orders of magnitude in mass 
(|Cappellari et all 1201 1[ ). In particular, 259 of the 260 
sample galaxies were searched for CO with the Insti- 
tut de Radioastronomie M illimetrique (IRAM) 30m tele- 
scope (|Young et al.l[20TTI ). The Combined Array for Re- 
search in Millimeter Astronomy (CARMA) interferome- 
ter is now mapping the brightest two thirds of the de- 
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tections in the J=l-0 transition of CO. Here we report 
on the SO field galaxy NGC 1266, the brightest of the 
CO detections in the ATLAS'^^ sample. The distance to 
NGC 1266 is taken from ATLAS^^, 29.9 Mpc, for which 
1" = 145 pc. The single-dish line profile requires a dou- 
ble Gaussian fit, with broad wings that exceed the escape 
velocity of NGC 1266. The interferometric CO observa- 
tions suggest that NGC 1266 possesses an extraordinar- 
ily dense, centrally compact molecular nucleus with a 
powerful outflow. In Sj2l we describe the molecular line 
observations using the IRAM 30m single-dish, CARMA 
and Submillimeter Array (SMA) telescopes. We show in 
§3 the two-component structure of the molecular gas - 
a massive central compact nucleus and a more extended 
outflow that appears to be escaping the galaxy, and de- 
rive the basic physical properties of the molecular gas. 
In §4, we compare NGC 1266 to other galaxies and dis- 
cuss the remarkable properties of its molecular gas. We 
summarize our main conclusions in §5. 

2. OBSERVATIONS 

NGC 1266 was observed as part of a flux-limited 
CO(l-O) and C0(2-l) survey of ah ATL AS^" galaxies 
carried out at the IRAM 30m telescope (jYoung et al.l 
120 lit ). The data consist of a single pointing at the galaxy 
center, covering a velocity width of 1300 km with a 
spectral resolution of 2.6 km s~^ at CO(l-O), centered 
on the systemic velocity of 2155 km s~^ a s listed in the 
LEDA catalog (based on absorption lines; iPaturel et al.l 
I2003| ). with a beamsize of 21'.'6 in the CO(l-O) line. The 
C0(2-l) covered velocity width of 1300 km s~^ with a 
spectral resolution of 5.2 km s^^, with a beamsize of 
12". A zcroth order baseline was subtracted, using line- 
free channels to estimate the continuum contribution to 
the flux. NGC 1266 has a peak main beam temperature 
Trr^h of 250 mK for CO(l-O) and 630 mK for C0(2-l) 
(lYoung et al. l [20Tll) . or 1.2 Jy for CO(l-O) and 3.5 Jy 
for C0(2-l) using a conversion factor of 4.73 Jy K~^, 
standard for the IRAM 30m. 

NGC 1266 was observed in CO(l-O) with CARMA 
between 2008 March and 2009 June in three different 
configurations: D-array (5" resolution at CO(l-O) with 
robust=0 weighting), B-array (0'.'7), and A-array (0'.'3). 
The primary beam had a diameter of 2' at CO(l-O), 
covering the full extent of the optical emission. Due to 
limited correlator bandwidth, a velocity extent of only 
420 km s~^ at 2.5 km s~^ resolution was available for 
the CO(l-O) line at the time of the observations, enough 
to adequately measure the properties of the center of the 
line, but not of the wings. D-array observations of the 
bright calibrator 0423—013 were taken in order to correct 
the slope across the bandpass, followed by alternating 
integrations on the source and a gain calibrator (0339- 
107), used to correct atmospheric phase fluctuations. 
The CARMA A and B arrays utili ze the Paired An - 
tenna Calibration System (FAGS; see lPerez et al.|[2010[ ). 
in which the 3.5 m dishes of the Sunyaev-Zel'dovich Ar- 
ray (SZA) are paired with CARMA antennas forming 
the longest baselines to correct for atmospheric phase 
fluctuations. The SZA antennas continually observed a 
bright quasar (in our case 0339—017) located 5?9 from 
the source, which was also used to calibrate the gains. 
The flux calibrator was 0423-013. At 29.9 Mpc, the 



highest angular resolution (0'.'3) corresponds to 44 pc. 
The data were reduced in the standard manner, using 
the Multichannel Image Reconstruction, I mage Analy- 
sis, a nd Display (MIRIAD) software package (jSault et al.l 
[1991 . 

SMA observations were also obtained between 2008 
September and 2009 March in two different configura- 
tions: compact (4" resolution at C0(2~l) with robust=0 
weighting) and extended (2"). The molecular species ob- 
served were C0(2-l) and i3cO(2--l) in the 230 GHz 
band, and CO (3-2) and HCO+(4~3) in the 345 GHz 
band. The primary beam had a diameter of 45" at 
C0(2-l) and 30" at CO(3-2), which includes ah emis- 
sion captured by the single dish. NGC 1266 was imaged 
with a bandwidth of 2600 km s~^ and a spectral reso- 
lution of 2.5 km s~^ at C0(2-l), centered on the sys- 
temic velocity, and a bandwidth of 1700 km s^^ and a 
spectral resolution of 0.7 km at CO(3-2), also cen- 
tered on the systemic velocity. We used J0423— 013 and 
J0238-f 166 as gain calibrators. The bandpass calibrators 
were 3C454.3 and 3C273, and the flux was calibrated 
with Titan and Ganymede. The zcroth order continuum 
emission of 40 mJy and 13 mJy at CO(3-2) and C0(2-l), 
respectively, was estimated using line-free channels and 
subtracted from the datacubes. These continuum levels 
are consistent with being due to the Rayleigh- Jeans tail 
of dust continuum. The data were reduced using the MIR 
package, and analyzed using MIRIAD. 

NGC 1266 was observed with the Expanded Very Large 
Array (EVLA) in the D configuration on 2010 March 20 
and 28, giving a total of 264 minutes on source with a 
correlator setup with 256 channels, each of 31.25 kHz 
width, centered on H I at the systemic velocity of the 
galaxy. The gain calibrator was J0323+0534 and the 
flux and bandpass calibrator was J0137-I-3309. Stan- 
dard observing and data reduction techniques were used, 
performed in the Astronomical Image Processing Sys- 
tem (AIPS) package. Observations were made at a fixed 
sky frequency, and the changing topocentric velocity cor- 
rection between the two observing dates amounts to 0.4 
channels; the AIPS task CVEL was used to shift one of 
the datasets before combining them. The data were Ban- 
ning smoothed in velocity, which helps to mitigate the ef- 
fects of interference at 1408 MHz. The final, continuum- 
subtracted image cube was made with robust weighting 
(robust parameter = 0), has a beam size 66'.'9x43'.'0, and 
covers the velocity range 1321 to 3014 km s~^ at a veloc- 
ity resolution of 13.4 km s~^. The rms noise level is 0.8 
mJy beam~^ per channel. There is no B I emission ev- 
ident anywhere within the primary beam (31' FWBM), 
but prominent B I absorption against the central contin- 
uum source in NGC 1266 is observed. 

Archival 1.4 and 5 GBz continuum data obtained at 
the Very Large Array (VLA) are available for NGC 1266 
at resolutions of 1" and 0'.'4, respectively. These data for 
project AB660 were originally taken in December 1992 in 
the VLA A configuration under the standard 2 x 50 MBz 
bands with left and right circular polarizations. Time on 
source in each frequency was 380 seconds. 3C84 was used 
as an absolute flux calibrator and B0 336-019 was the 
gain calibrator (|Baan k K16cknei|[2006l ). The data were 
reduced in the standard manner using the AIPS package. 

3. RESULTS 
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Figure 1. (Left:) CO(l-O) integrated fluxes from tlie IRAM 30m (black and gray) and CARMA A, B & D-arrays (red). Tiie total 
CARMA A, B & D velocity coverage at CO(l-O) is 390 km s^^, compared to 1300 km s^^ for the IRAM 30m. The total velocity widths of 
both observations are present. Dashed red lines indicate the CARMA bandwidth. (Middle:) CO(2-l) fluxes from the IRAM 30m (black) 
and the combined SMA compact (C) and extended (EX) conflgurations (blue). (Right:) CO(3-2) flux from the combined SMA C and EX 
configurations. CARMA recovers 20% more flux than the single-dish, and the SMA 30%, within single-dish baseline and flux calibration 
uncertainties. 



3.1. Molecular Line Profile 

Figure [T] presents the integrated spectra of the ^^CO 
data used in this paper, showing good agreement be- 
tween the CO(l-O) and C0(2~l) fluxes recovered from 
the single-dish and the interferometers. Also of note is 
that the CO (3-2) line profile from the SMA has the 
same shape as the C0(2-l), including wing emission. 
The CARMA observations recover « 120% of the single- 
dish flux, and the SMA observations w 130%. Errors in 
the single-dish baselines and the intrinsic flux variations 
of the mm-calibrators can account for the discrepancy. 
Therefore, it appears that both CARMA and the SMA 
recover all of the single-dish flux. 

With its narrow single-peaked central velocity compo- 
nent (hereafter CVC) and high-velocity wings, the in- 
tegrated profile shape is mor e akin to t hose observed 
in protostellar outflows (e.g. iBallv &: Lad a 1983) than 
those of typical external galaxies, where a double-horned 
profile, is usually seen. 

Fig. [2] shows the fits and residuals of the single and 
nested Gaussians. Table 1 gives a summary of the fitted 
Gaussian parameters. We find that the broad wing and 
CVC, which we identify respectively with the broad and 
narrow Gaussian of the nested Gaussian fits, encompass 
respectively 34% and 66% of the C0(2-l) single-dish in- 
tegrated flux, which has the best signal-to-noise ratio. 
When comparing the line to the escape velocity radial 
profile measured in the plane of the galaxy (see Figure 
[S]), we find that at least 2.5% of the wing emission is from 
velocities exceeding the innermost local escape velocity 
(calculated at 2") of Vcsc ~ 340 km s~^. The true per- 
centage may be higher, depending on the geometry of the 
high velocity gas. If the molecular gas is located either 
out of the plane of the galaxy or at larger radii, the local 
potential (and thus Vcsc) would be much reduced. The 
Vcsc profile was calculated using Jeans Anisotropic Multi- 
Gaussian Expansion (MCE) modchng (JAM; ^33)) . to 
extract the true potential of the galaxy by mo deling the 
observed stellar kinematics (jScott et al.ll2009l in prep). 



Table 1 

Fitted Gaussian parameters for the CO(l-O), CO(2-l) and 
CO (3-2) lines 





CVC 


broad wing 


CO(l-O) 








2162 ± 1.4 km s"! 


2142 ± 20 km s"! 


FWHM = 


114 ± 3.9 km s-i 


353 ± 17 km s"! 


rmb(peak) = 


0.21 ± 0.006 K 


0.03 ±0.005 K 


Icq = 


25.5 ± 1.1 K km s^^ 


9.6 ± 1.8 K km s^^ 


CO(2-l) 








2161 ±0.5 km s"! 


2151 ±4.4 km s^i 


FWHM = 


128 ± 1.9 km 


353 ± 17 km s"! 


rmblpeak) = 


0.53 ± 0.01 K 


0.10 ±0.01 K 


Icq = 


71.1 ± 1.6 K km s-i 


36.7 ±3.8 K km s-^ 


CO(3-2) 








2166 ±0.4 km s"! 


2162 ± 3.7 km s"! 


FWHM = 


134 ± 1.3 km s-i 


353 ± 17 km s"! 


^peak 


4.15 ± 0.03 Jy 


0.82 ±0.03 Jy 


Sco = 


589 ± 7.5 Jy km s^l 


306 ± 18 Jy km s^^ 



Note: Ico is the integrated CO line intensity from the IRAM 
30m, similarly for Sco a^nd the SMA. The broad wing compo- 
nent accounts for 27% and 34% of the total flux for CO(l-O) and 
CO (2-1), respectively. The broad wing component of CO (3— 2) also 
contributes 34% of the total flux. The 7% difference in the broad 
wing contribution between the transitions is attributed to the un- 
certainties in the fits to the line profiles. We will adopt 34% as 
the more likely value, due to the better signal-to-noise ratios at 
CO (2-1) and CO (3-2). CO (3-2) data were only obtained with an 
interferometer, thus the natural unit of this line Jy. To convert 
the CO (3-2) fiuxes to Kelvin, we first divide by 9 (the ratio of the 
CO(3— 2) and CO(l-O) frequencies squared), then multiply by the 
K per Jy conversion factor of the adopted beam of 4.73 Jy K~^. 
The errors listed are the formal uncertainties of the fits. The value 
and uncertainty of the FWHM of the broad wing component for 
all transitions are adopted from the fit to the CO(2-l) line. 

3.2. Molecular Gas Mass 

Using the CO (1-0) single-dish flux quoted in Young 
et al. (2011), a column density of molecular hydrogen 
N(H2) = Xco J Tmh dw, where the conversion factor 
Xco 2 X 10^" cm-2 (K km s'^)''^ (iDame et al.|[200l . 
a beamsize of 21'.'6, and a beam area of 529 sq.arcsec 
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Figure 2. CO(l-O) (top), CO(2-l) (middle) and CO(3-2) (bot- 
tom) integrated profiles of NGC 1266 from the IRAM 30m and 
SMA (CO(3-2) only). We also show a comparison between a one 
Gaussian fit the the data (red) and a two-Gaussian fit (blue) con- 
sisting of a narrow (blue dot-dashed line) and a broad (blue dashed 
line) component. Residuals are plotted below the fits and show that 
a two-Gaussian fit is a much better match to the data. The width 
of the broad component is fixed to 353 km s~^, initially constrained 
using the fit the to the CO(2-l) line. The goodness of fit improves 
significantly with the addition of the broad component, illustrated 
both in the measured ^ well as in the residuals. 



(whic h includes a bcamshapc correction of 1 /ln(2) ; iBaarsI 
l2007f ). the total mass of the molecular gas would be 
1.7 X 10^ Mq, including a correction factor of 1.36 for 
He. However, the wing emission is optically thin (see 
i i3.4[) and contributes little to the total mass. Consider- 
ing only the CVC, the molecular gas mass is 1.1 x 10^ 
Mp), including He. 

iDraine et al.l ()2007D derive a total dust mass of 1.02 
X 10'' Mq in NGC 1266 (rescaled to our distance of 
29.9 Mpc), using Spitzer Infrared Nearby Galaxy Sur- 
vey (SINGS) data. Since a significant H I reservoir 
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Figure 3. The radial escape velocity profile of NGC 1266. The es- 
cape velocity was derived using Jeans Anisotropic Mu lti-Gaussian 
Expa nsion modeling (JAM) of the stellar kinematics IjScott et al.l 
[20091 . 2010, in prep). 

is not observed in NGC 1266 (see below), we may as- 
sume that all of the measured dust is associated with the 
molecular gas. Adoptin g a gas-to-dust mass ratio of 100 
(jSavage fc Mathidll979[ ). the dust mass implies a molec- 
ular gas mass of 1.0 x 10^ Mq. The derived dust mass 
is thus in good agreement with the mass derived from 
the CO measurements, and a standard Xqo conversion 
factor. 

3.3. Core of the Emission 

Figure |4] shows the full spatial extent of the CARMA 
CO(1~0) integrated emission {Av = 390 km s~^), over- 
laid on an R-band image of the galaxy from SINGS, il- 
lustrating that the vast majority of the molecular gas is 
restricted to the innermost region of NGC 1266. The 
righthand panel of the figure shows that the molecular 
gas is not well described by a single spatial component, 
but instead as both a diffuse, relatively extended com- 
ponent, which we will call the "envelope," and a highly 
concentrated nuclear structure, which we will call the 
"nucleus" . 

In order to separate the nucleus from the envelope in 
the CARMA map, we take a slice across the integrated 
intensity map at a position angle of 145°, to trace the 
elongation of the narrow component, and fit two Gaus- 
sianS; to the nucleus and the envelope, respectively. Fig- 
ure [5] illustrates the radial cut as well as the the nuclear 
Gaussian component, which has a characteristic FWHM 
of 1'.'05. We use FWHM/2 to define a characteristic de- 
convolved radius i?nuc of « 60 pc (after accounting for 
the CARMA beam). Using this, the peak intensity of 
the profile in Fig. [5] (24.3 Jy beam^^ km s~^), the pixel 
size of 0"08 pix-^, and Xqo = 2 x 10^° cm^^ (K km 
s~^)~^, we derive a molecular gas mass for the nucleus 
of Mnuc,H2-i-Ho = 4.1 X 10^ M© , including He. As the spa- 
tial distribution of the broad wing emission in the central 
region of the CARMA CO(l-O) integrated image is un- 
known, we do not subtract a potential wing contribution 
from this nucleus molecular gas mass. If we assume that 
the rest of the mass in the CVC comes from the envelope. 
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Figure 4. (Left) Contours of the integrated CO(l-O) molecular gas detected by CARMA (blue), overlaid with the isophotes (black) and 
grayscale of an R-band image from the Cerro Tololo Inter-american Observatory (CTIO) 1.4m telescope (SINGS). (Right) A zoomed-in 
version of the CO(l-O) integrated intensity map. The synthesized beam (0"7x0'.'6) is plotted in the lower right corner. The CARMA 
CO(l-O) integrated intensity map was created by summing the full velocity width of the CARMA cube, including the A, B and D arrays. 
Robust = weighting was used in order to reveal low surface brightness, more diffuse structures at the expense of resolving the most 
compact components. Contour levels are 3, 6, 9, 12, 15, 18, 21, 24 and 27 Jy beam^^ km s^^ (rms = 1.0 Jy beam^^ km s^^). 



we deduce that Monv ~ 7.0 x 10^ M©, including He. 
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Figure 5. Slice of the CO(l-O) integrated intensity map from 
CARMA at a 145° position angle (see inset). The radial profile re- 
quires two Gaussians for a good fit, an "envelope" and a "nuclear" 
component. The nuclear component (red dashed line) has a peak 
intensity of 24.3 Jy beam^i km s~i and a FWHM of 1.05", (0'.'82 
after correcting for convolution effects from the beam). The radial 
profile of the beam is also shown (red dotted line), to illustrate 
that the nucleus is spatially resolved. 

Figure [6] is a map of the first moment of the CO(1~0) 
emission, and shows a clear velocity gradient across the 
nucleus. Figure [7] is the position- velocity (PV) diagram 
at a position angle of 90° (i.e., roughly along the kine- 
matic major axis), suggestive of a rotating disk with a 
velocity gradient of 1.8 km s~^ pc~^. Rotation can be 
traced on the PV diagram to a radius of 54 pc, consistent 
with the radius of the nucleus derived above. At i?nuc ~ 
60 pc, we derive a maximum projected rotation velocity 



of the nucleus of Urot ~110 km s ^. The corresponding 
enclosed dynamical mass, assuming circular rotation and 
spherical symmetry (Afdyn,nuc = i?nucWrot/ sin^ inucG, 
where inuc is the inclination of the nuclear disk and 
G is the gravitational constant), is uncertain because 
of the unknown inclination, but the formal lower limit 
(ijj^^ = 90°) is 1.5 X 10® M0. A rough inchnation esti- 
mate based on the axial ratio of the CO(l-O) integrated 
intensity map shown in Figure |4] and the assumption of a 
thin disk yields inuc ~ 34°, and thus an enclosed dynami- 
cal mass Mdyn.nuc ~ 4.9 x 10® Mq. The dynamical mass 
calculated using the rotation of the nucleus is therefore 
consistent with the ^^CO-derivcd molecular gas mass of 
the nucleus, when taking into account a correction asso- 
ciated with the inclination. The stellar contribution to 
the nucleus mass is discussed later. 

An independent estimate of the mass of molecular gas 
in the core was obtained via dynamical modeling of the 
stellar kinematics. For this we constructed a MGE model 
lEmsellem et all ()1994D of the stellar mass distribution of 
NGC 1266 from our own R-band Isaac Newton Telescope 
(INT) pho tometry (Scott et al. 2011, in prep) and the 
software of lCappellanl (|2002D . We used this MGE model 
to build an axisymmetric dynamical model of the stel- 
lar ki nematics based on the JAM formalism ([CapDellaril 
|2008[ ). The global anisotropy, inclination and mass-to- 
light ratio (M/L) of the galaxy were fitted to the large- 
scale ATLAS^° SAURON stellar kinematics (CappcUari et 
al. 2011b, in prep) and will be discussed elsewhere. In 
this particular case the JAM model also explicitly in- 
cludes as an extra free parameter: the mass of the molec- 
ular gas in the core, modeled as a dark Gaussian with a 
width of 1'.'2. The mass of molecular gas that best fits the 
SAURON stellar kinematics in the centre is Mgas « 8 x 10® 
Mq, with a factor « 3 uncertainty dominated by sys- 
tematic effects. This value is an order of magnitude 
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larger than the mass in stars within 100 pc of the center, 
-^*,ioopc ~ 5 X 10^ Mq, as determined from an analo- 
gous MGE mass model without the molecular core. This 
result is thus consistent with the mass in the center of 
NGC 1266 being dominated by molecular gas. 




3''16'"00.8^ 00.7' 
RA (J2000) 



Figure 6. CO(l-O) mean velocity map from CARMA overlaid 
with contours from the corresponding integrated intensity map 
(black and white). A velocity gradient is clearly seen across the 
nucleus, interpreted as a thick rotating disk (see text). These ve- 
locity moments were constructed from the CARMA A and B array 
data only (synthesized beam 0'.'5x0'.'4), to extract structures at 
the smallest spatial scales at the expense of the extended, diffuse 
emission. The velocity contours are spaced at 25km s~^, or 2069, 
2094, 2119, 2144, 2169, and 2194 km s'^. The systemic velocity 
of 2160 km s~^ is denoted with an arrow. T he centroiding accu - 
racy is derived using the FWHMboam/2 SNR UTavlor et al.lll999l) . 
where SNR is the signal-to-noise ratio and is at least 3 per channel, 
yielding an accuracy of at least O^'OS. 

One final consistency check on the molecular gas 
mass can be performed using the CO (1-0) single- 
dish profile obtained as part of the dense gas census 
of ATLAS^° galaxies (Crocker et al. 2011, in prep), 
shown in Fig. M The total i3cO(l~0) line inten- 
sity is 0.98 K km s"^ within the 22'.'4 beam of the 
IRAM 30m telescope. We must first correct /(^^CO) 
to reflect a more realistic size for the emitting region, 
and adopt a radius of ~ 1". The correction corre- 
sponds to the ratio of the two solid angles subtended: 
^hcam/T^R^, in this case 567 arcsec^/3.14 arcsec^. The 
resulting corrected I{^'^CO) is 177 K km s~^. Using 
the ^'^CO to Ay conversion from iLada et al.l (fTOOl . 
/("CO)/(K km s'l) = 1.88 + 0.72 mag-^ Ay results in 
an extinction Ay ~ 243 mag. Converting Ay to 
N(H2) using N (li2 ) / Ay = 9 x 10^° cm"^ . i ^^g"^ from 
ISchultz fc Wie^ mm) and IBohlin eFall (flOTl . and 
then converting N(H2) to a molecular gas mass, yields 
M„uc « 3.2 X 10® Mq, corrected for He. Therefore the 
^■^CO-derived mass is in reasonable agreement with the 
^^CO-derived molecular gas mass. 

With a molecular gas mass of 4.1 x 10® M© (including 
He), assuming a circular disk, the mean molecular hydro- 
gen surface density in the nucleus is ~ 2.7 x 10* M© 
pc~^ (corresponding to a mean column density iV(H2) ~ 
1.7 X lO^'' cm~^) and the mean volume density of H2 is 
at least 6.9 x 10'^ cm~'^ (lower limit assuming spherical 




Offset (arcsec), p.a. = 90' 

Figure 7. Position-velocity diagram of NGC 1266 for our 
CO (1-0) A and B array CARMA data, taken at position angle 
of 90° . The observed lack of a turnover is indicative of the gas not 
reaching the flat part of the rotation curve. The data are consis- 
tent with rotation with a gradient of 1.8 km s~^ pc~^ in the inner 
parts. The red line is a trace of the peak at each velocity. The 
dark red dashed lines indicate v^ys and the center. 

symmetry). Given that tracers with higher critical den- 
sities, such as i^CO, CS, HCO+ and HCN, are robustly 
detected in the nucleus (Alatalo et al. 2011, in prep), the 
molecular gas is likely to be clumpy. 

The dynamical mass, JAM mass and ^'^CO-derived 
mass agree quite well with the ^^CO(l-0)-derived molec- 
ular gas mass of the nucleus, indicating that using the 
standard Milky Way conversion factor Xqo = 2 x 
10^° cm~^ (K km s"-'^)"^ is appropriate to determine 
the mass of the molecular material in the center of 
NGC 1266. The JAM modeling also indicates that the 
stellar contribution to the mass of the region with molec- 
ular gas is negligible. 

3.4. Wings of the Emission 

We use the C0(2"l) SMA data to resolve the spatial 
extent of the high velocity wings of the molecular gas, 
because of the limited velocity coverage of the correlator 
when the CARMA data were taken. In order to minimize 
contamination from the nucleus, we made a map of the 
emission over the velocity ranges v=1590-1890 km s~^ 
and u=2350"2650 km s~^, eliminating the central 460 
km s~^. Figure [9^ shows the velocity ranges of the CVC 
and the wings of the line used for mapping. The offsets 
from Vgys were chosen to achieve equal signal-to-noise ra- 
tios of 6 for the peaks of the red- and blue-shifted images 
in Figureint, which shows that the wings are considerably 
more extended spatially than the nucleus. If we model 
the outflow as a b i-spherical structur e, similar to what is 
observed in M82 (jWalter et al.ll2002D . and represent the 
red- and blueshifted lobes by two tangent spheres whose 
point of contact is the nucleus, the distance the molecu- 
lar gas has traveled from the center is equal to the width 
of the lobes. We measure the width of the blueshifted 
and redshifted lobes to be « 3'.'4 and 4'.' 2, respectively. 
Taking the average of the lobe extents, and correcting for 



AGN-drivcn Molecular Outflow from NGC 1266 



7 



convolution with the 2" SMA beam gives us a character- 
istic extent of w 3'.'2, or 460 pc. The red- and blue-shifted 
wings are also spatially offset from one another. The po- 
sition centroids of the lobes are separated by w 2'.'3, thus 
each is offset by about 170 pc from the galaxy center. 
The large spatial extent of the wing emission rules out 
an interpretation whereby the high gas velocities are the 
result of gas having fallen into the central potential well 
of the galaxy. We thus conclude that the source of the 
high velocity wing emission is a kpc-scale molecular out- 
flow, with i?outfiow ~ 460 pc. 

Further evidence that the wings represent an outflow 
is that the axis connecting the lobes is intersected by the 
CVC and is not aligned with the kinematic major-axis of 
the nuclear gas disk. This suggests that the broad wing 
emission is composed of gas that is being expelled, and 
is not an extension of the nuclear material. If we esti- 
mate the inclination of the outflow (with respect to the 
plane of the sky) using the average offset of the centroids 
with respect to the nucleus (2'.'3) divided by the average 
extent of the lobes (3'.'2) calculated above, we obtain an 
inclination angle of roughly 20°. 

The mass of outflowing gas is somewhat uncertain 
because the wings are unambiguously detected only in 
^^CO (Fig. [T|), requiring that the density, optical depth 
and temperature of the outflow be determined from the 
excitation of the three lowest CO rotational transitions. 
We use the beam-corrected (to the single-dish CO(l-O) 
beam of 21'.'6) intensities of the wing emission from the 
CO(l-O) and (2-1) IRAM 30m single-dish spectra and 
the CO(3-2) SMA spectrum (see Table 1). The beam- 
corrected values are 11.93 K km s~^ for CO(l-O), 9.16 
K km s-^ for C0(2~l) and 7.16 K km s'^ for CO(3-2). 

We use the RADEX large vel ocity gradient (LVG) soft- 
ware (jvan der Tak et"aLl2007t ) to determine volume den- 



sities and column densities in the wind. With RADEX, 
assuming a line width Wfwhm = 353 km s~^, we find that 
the conditions required to reproduce the beam-corrected 
flux ratios above are an H2 density n(H2) « 10^ cm~^, 
a kinetic temperature Tkin ~ 100 K, and a CO col- 
umn density A^(CO) 1.0 x 10^^ cm'^. The iV(CO) 
w as also constrained usi ng on the optically thin estimate 
of IKnapp fc Jural (119761) assuming subthermally excited 
CO, and the result is consistent with that from RADEX. 
Assuming a CO/H2 abundance ratio of ID""', the derived 
A^(CO) and our adopted beamsize of 21'.'6 imply a molec- 
ular outflow mass of 2.4 x 10^ Mq, including He. It is 
however possible that the outflowing molecular mass is 
higher. The optically thin approximation (i.e. that we 
are seeing all CO molecules) provides us with a lower 
limit to the total outflowing molecular mass. If the state 
of the molecular gas in the outflow of NGC 1266 is more 
akin to gas in either ULIRGs or giant molecular clouds 
(GMCs) in the Milky Way, where a conversion factor 
(Xco) is used, the total outflowing mass could be a fac- 
tor of w 10 - 20 higher. 

3.5. Additional Data 

While no H I emission is detected in NGC 1266, H I 
appears in absorption against the 1.4 GHz continuum 
source. The H I observations have thus far only been 
taken at low resolution, so are presently unresolved. The 
top panel of Figure [TOl shows the H I absorption profile, 
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Figure 8. ^^CO(l-O) spectrum taken with tiic IRAM 30m tele- 
scope, as part of a larger census of dense gas in ATLAS'^^ ETGs 
(Crocker ct al. 2011, in prep). The beam at this frequency has a 
22" half-power beam width (HPBW) and the channels shown are 
30 km s-i. 

which exhibits a broad blueshifted velocity component, 
and the fitted Gaussians describing the absorption pro- 
file. The bottom panel of Figure ITUl shows a good corre- 
spondence between the blueshifted absorption feature in 
the H I and the high-velocity wing of the C0(2-l) emis- 
sion line profile. This identifies the blueshifted molecular 
gas as moving out (as opposed to infalling from behind), 
and implies that the outflow in NGC 1266 has multiple 
phases. 

Calculating the column density from the absorption 
feature, assuming Tspin = 100 K, reveals that the H I 
column in front of the continuum source is A^h = 2.1 x 
10^^ cm~^, with the outflowing component contributing 
8.9 X 10^° cm~^ and the systemic component contributing 
^. If we assume that the H I and CO are 
^outflow = 460 pc), and that the total 



1.2 X 1021 cm 
co-spatial (i?Hi 

mass of outflowing H I is twice what we calculate (since 
we are only able to detect the blueshifted lobe) , we derive 
the total H I mass in the outflow to be 4.8 x 10^ Mg 
per lobe, totaling 9.5 x 10^ Mq. If we include the H I 
contribution, the total neutral gas mass (H I -I- H2 + He) 
of the outflow is thus 3.3 x 10^ Mg. 

The molecular gas (as traced by CO) and atomic gas 
are likely not the only constituents of the mass in the 
outflow. By not accounting for other states of the gas, 
we are underestimati ng the true outflow mass. For ex- 
ample, iRoussel et al.l 



(|2007[ ) detect a large reservoir of 
warm H2 in NGC 1266 using Spitzer, with a total mass 
of il'/Hj.warm ~ 1-3 X 10^ Mg . It is vcry likely that 
this mass of warm H2 belongs to the outflow rather 
than the nuclear disk. The total warm H2 luminosity 
is iHa.warm « 1 X CTg . FoT the uuclcar region 
to sustain this luminosity of warm H2 would require an 
AGN with a bolometric luminosity of > 2 x 1 0'^'^ erg s~^. 
assum ing flrst that Lx agn ^. lOO^q -^ warm, lOgle et ahl 
[2OT0t and Lboi,AGN > 16£x, E^EQQ^,' a factor of three 
large r than Lfir (7 x 10^"^ erg s~^; iGil de Paz et al.l 
[200% . If this warm H2 is indeed part of the outflow, 
our estimate of the total mass of the outflow should in- 
crease by 40%. Adding in the contribution of the ionized 
and hot gas would increase the mass of the outflow fur- 
ther, thus we consider our estimate of 3.3 x 10'' Mg to 
be a conservative lower limit. 
Combined with the molecular emission and H I absorp- 
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Figure 9. (Top:) CO(2-l) spectrum from the IRAM 30m tele- 
scope, indicating the velocity ranges used. (Bottom:) CO core 
and wings in NGC 1266, overlaid on a greyscale Ha narrow-band 
image from SINGS. Superimposed contours are from the CARMA 
CO(l-O) integrated intensity map (yellow) and the SMA CO(2-l) 
redshifted (red) and blueshifted (blue) wings. Contours are 3, 4, 5 
and 6 Jy beam~^ km for the red and blue components (rms 
= 1 Jy beam-l km s"!) and 6.4, 8.6, 10.8, 13.0, 15.2 and 17.4 Jy 
beam~^ km s^^ for the yellow component (rms = 1.1 Jy beam~^ 
km s"-*^) 



tion data, X-ray, Ha, and radio continuum data further 
elucidate the nature of the outflow. Figure [Tl] shows 
a comparison o f the Hct narrow-band image from the 
SINGS survey (iKennicutt at all 12001 a 1.4 GHz ra- 
dio continuum map (jBaan fc Klockneii [20061 ) . and the 
unsmoothed X-ray map from the Chandra X-ray Tele- 
scope Advanced CCD Imaging Spectrometer (ACIS) on 
the same spatial scale (Alatalo et al. 2011, in prep). A 
spectral fit to the Chandra data co-added over the extent 
of the emission shows that the X-ray emission is domi- 
nated by thermal bremsstrahlung. The 1.4 GHz con- 
tinuum shows both a centrally peaked, unresolved core 
(possibly from a radio jet), as well as two more diffuse 



spurs running from the southeast to the northwest. The 
lobed structure first apparent in the molecular gas (Fig. 
[9]) is seen in all these tracers. In fact, the X-rays and the 
1.4 GHz spurs, which do not suffer significantly from ex- 
tinction, show an asymmetry between the foreground and 
background lobes, relative to the CO at Wgys- This likely 
means that the relative faintness of the Ha flux from the 
redshifted lobe as compared to the blueshifted lobe is not 
primarily due to extinction. The spatial relationship be- 
tween the molecular gas and the X-ray, 1 .4 GHz and Ha 
emitting material supports this picture. The Ha, radio 
spurs, and X-rays likely are emitted from the interface 
where the outflowing material collides with the interstel- 
lar medium (ISM) of NGC 1266. This interpretation 
is supported by the thermal bremsstrahlung emission in 
the Chandra spectrum, although the radio spur emission 
could also possibly come from a radio jet. 

The geometry of the CO outflow can be explained us- 
ing two possible scenarios. Firstly, it could be due to 
a molecular outflow launched deep in the nucleus of the 
galaxy, confined azimuthally by the massive molecular 
disk. Secondly, it could directly originate from the molec- 
ular disk, and be entrained by hot winds launched by the 
AGN, surviving only because of the extreme density of 
the molecular disk. In either of these scenarios, at least 
some of the molecular material escapes from the galaxy 
still intact, and will energize the intergalactic medium 
(IGM) . This is especially true if the outflow continues to 
accelerate from its launch point. 

3.6. Derived Properties 

A summary of the derived properties of NGC 1266 is 
presented in Tabled The mean molecular hydrogen sur- 
face density of the nucleus is 2.7 x 10^ pc~^, corre- 
sponding to a mean column density Af(H2) ~ 1.7 x 10^^ 
cm~^, two orders of magnitude higher than the surface 
density of individual GMCs in the Milky Way, and com- 
parable to what is seen in some of the most luminous 
starburst galaxies (Kennicutt 1998). 

For an outflow with a mass of 3.3 x 10^ Mq and an out- 
flow velocity of half the full width-half maximum (ufwhm) 
of the broad wing component (177 km s^^), the kinetic 
energy associated with the outflow is 1.0 x 10^^ ergs, 
equivalent to the total kinetic energy expelled from 10^ 
supernova explosions. The dynamical time, Tdyn, for the 
outflow to reach its maximum radius of 460 pc at this 
velocity is only 2.6 Myr. If the edge of the CO emission 
in the lobes is produced from only the highest velocity 
gas, the relevant timescale is even shorter, Tdyn 1 Myr. 

From the total kinetic energy and Tdyn, the mechanical 
luminosity of the outflow Loutfiow = 1.3 x 10*^ erg s~^. 
The mass outflow rate is M = Moutfiow/Tdyn = 13 Mq 
yr"^. At this rate, if the source of the gas in the outflow 
is the nuclear molecular disk, the nucleus of NGC 1266 
will be completely depleted in Tdop = 32 Myr, or 85 Myr 
for all of the gas in the center (nucleus -|- envelope) , both 
very short timescales. 

4. DISCUSSION 

4.1. Uniqueness of NCC 1266 

NGC 1266 is unusual in several respects. There is 
a large reservoir of molecular gas concentrated within 
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Figure 10. (Top:) Continuum-subtracted H I profile in the 
central pixel of the unresolved map observed with the EVLA D- 
array (black). The continuum level for the source is 113 mjy. 
The blueshifted wing is evident in the absorption profile, reaching 
> 200 km s~^ before dropping into the noise. Overlaid (red) is a 
two-Gaussian fit to the absorption profile, including a blueshifted 
high-velocity contribution (dashed) and a deeper absorption con- 
tribution near Vsys (dot-dashed). (Bottom:) A direct comparison 
between the H I absorption profile (black) and the CO (2-1) emis- 
sion spectrum from the IRAM 30m telescope (blue dotted). The 
CO emission and H I absorption profiles trace each other well at 
high blueshifted velocities, and the marked absence of redshifted 
H I absorption is confirmation that NGC 1266 harbors an H I 
outflow. 

the central 100 pc. The concentration of H2 in the nu- 
cleus implies that the gas must have lost nearly all of 



Table 2 

NGC1266 Characteristic Properties 



-^nucleus 


4.1 X 10** Mq 


-^nucleus 


60 pc 


Mcvc 


1.1 X 10^ Mq 


< S(H2)nuclcus > 


2.7 X 10* M0 pc-2 


< Af(H2)nuclcus > 


1.7 X 10^4 cm-2 


< 'l(H2)nuclcus > 


6.9 X 10^ cm"3 


I+H2, outflow 


3.3 X 10'^ M0 


.outflow 


2.4 X 10'' Mq 


-^outflow 


450 pc 


^outflow 


177 km s~^ 


Tdyn 


2.6 Myr 


M 


13 Mq yr-i 




1.0 X 10^^ ergs 


-^outflow 


1.3 X 10''! erg s"! 


^dcp 


85 Myr 



its angular momentum. Although such high central con- 
centrations of H2 are observed in interacting galaxies, 
NGC 1266 is an isolated galaxy with no sign of interac- 
tion or merger. 

Is the outflow in NGC 1266 a common evolutionary 
stage through which most ETGs pass, or is NGC 1266 
a pathological phenomenon? The short gas depletion 
time of w 85 Myr is reasonably consistent with finding 
only o ne such case so far in the 60 ATLAS'^'^ CO detec- 
tions (|Young et al.ll20l"TI ). The expectation would be to 
find 3, assuming that the typical molecular gas deple- 
tion time due to star formation in E TGs is the sa,me as 
the 2 Gyr found i n late- type galaxies (jLerov et al.ll2008l : 
iBigiel et al.ll2008D . However, the sheer strength of the 
CO line in NGC 1266 certainly aided the outflow detec- 
tion, and it is likely that, given the lower signal-to-noise 
ratios present in the rest of the sample, the low surface 
brightness features at high velocities would generally be 
lost in the noise. 

Optically, NGC 1266 does not stand out within the 
ATLAS'^'^ sam ple. It is a fast rotator with regular stel- 
lar kinematics (jKrainovic et al.l 120 lit ) , and is of normal 
metallicity. Its absolute magnitude, Mk = —22.9, cor- 
responds to a stellar mass just below the median of the 
sample. On the other hand, its radio fiux at 1.4 GHz 
places NGC 1266 in the top 10% of radio emitters within 
the sample, although this emission could be explained 
by the FIR-radio SF correlation based o n the amount 
of m olecular gas detected in the source (jMurgia et al.l 
120051). Among the - 100 known ETGs with detected 



molecular gas (iWiklind et al.lll995l: 


Welch & Sagd 120031: 


ISage et al.ll2007llCombes et al.ll2007 


; Young et aL 2011), 



w 2 have been previou sly mapped interferometrically 
(see lCrocker et aDl201lL and references therein). Com- 
pared to those, the molecular gas in NGC 1266 stands 
out both in terms of its compactness and its large mass. 
The star-forming ETGs seem to exhibit only a narrow 
range of mean molecular gas surface densities, 10-500 
Mq pc~^, similar to the range of Milky Way values. The 
surface density of the nucleus of NGC 1266 (« 2.7 x 10'' 
Mq pc~^) is two orders of magnitude larger than this. 

NGC 1266 should perhaps be compared to radio 
galaxies, 1/4 - 1/3 of which have been sh own to con- 
tain si gnificant re servoirs of molecular gas (jEvans et al.l 
120051 lOgle et all [2010). The H I absorption pro- 
file seen in NGC 126(3 is akin to that in the hand- 
ful of known ra dio galaxies with outflowing H I 
IC 5063. 3C 305 (iMOTganti et al.l lT998t lOosterloo et al 
2000 HMOTganti et al.ll2005allbll. 4C 12.50 (iMorganti et al 



200' 



and" 3C 293 (jMorganti et all l2003t lEmonts et al 
2001 . with mass outflow rates ~ 50 Mq yr ^ at high 
(> 1000 km s~^) velocities. The H I in each of these 
galaxies is detected in absorption against a radio source, 
and covers hundreds of parsecs to kiloparsecs. While the 
outflowing H I mass in NGC 1266 is orders of magnitude 
lower than in these galaxies, it is possible that both types 
of systems share a common origin. However, although 
the luminosity and structure of the H I absorption do 
mirror those in rad io ga laxies, the sample s described by 
lEvans et all (|2005f) and lOgle et all (|2010D and the out- 
flowing systems all show signs of interactions, compan- 
ions and tidal features, which are all markedly absent 
from NGC 1266. 
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Figure 11. (Left:) A RGB image of the Hq (red), radio continuum (green), and Chandra X-ray (blue) emission. A spatial correlation 
is clearly seen in the area of the blueshifted CO outflow, revealing that these features are co-spatial. We conclude that the Ho, the spurs 
in the radio continuum image, and some of the X-ray emission originate from the outflowing material, although the 1.4 GHz continuum 
emission could also be co ming from a radio jet. ( Top right:) Hq image from the SINGS survey. (Middle right:) VLA A-array 1.4 
GHz continuum emission I IBaan fc Klockrierll2006D . The unresolved peak in the emission is assumed to be from an AGN, and the lobe 
extending southeast to northwest is assumed to trace the interface between the outflowing material and the galaxy ISM. (Bottom right:) 
Un-smoothed Chandra X-ray image (Alatalo et al. 2011, in prep). The majority of the X-ray photons from NGC 1266 can be fit with a 
thermal bremsstrahlung spectrum with an excess of hard X-rays. All images are at the same spatial scale. 



4.2. Star Formation in NGC 1266 

Given the high surfaee density of the molecular gas, 
NGC 1266 ought to be forming st ars. Based o n 
the Kennicutt-Schmidt (K-S) relation (|Kennicuttlll998D . 

~ 2.7 X 10"^ Mq pc~^, and a radius of 60 pc, the 
expected star formation rate (SFR) of the nucleus alone 
is 3.1 Mq yr~^, adjusted for a Kroupa initial mass func- 
tion (IMF). In order to understand the star formation 
efficiency of the gas, we use multiple tracers to place con- 
straints on the SFR, including the total far-infrared lumi- 
nosity (Xfir), the 24/.tm emission, and the 8/^m emission 
from polycyclic aromatic hydrocarbons (PAHs). Those 
three tracers are systematically affected by the presence 
of an AGN, but in different ways, so we use all three to 
define a range of possible SFRs. Both the 24/xm emission 
and LpiR suffer from contamination by a buried AGN, 
and will thus systematically overestimate the SFR. PAH 
molecules are destroyed by the harsh ionization fields as- 
sociated with strong A GN, and wi ll thus possibly under- 
estimate the true SFR ()Voitlll992D . 

If we assume that all of the far-infrared (FIR) lu- 
minosity (7 x 10^ erg s"!; iGil de Paz eTall [20071: 
iCombes et al.l 120071 ) is from star formation (an up- 
per limit since the nucleus of NGC 1266 contains an 
AGN), we calculate a global obscured SFRfir of 2.2 
Mq yr~^ using iKennicut^ (I1998D . re-calibrated for a 
Kroupa IMF. The unobscured SFR from the measured 
far ultraviolet (FUV) contributes negligibly to the over- 



ah SFR (SFRfttv = -003 Mq yr'^; IGil de Paz eFall 
I2007t iLerov et al.il2008D . We also calculate the SFR us- 
ing the 24/xm and PAH emissio n, measured by Spitze r 
as part of the SINGs survey (jKennicutt et al.l |2003() . 
The SFR is related to the 24/xm luminosity (L24) fol- 
lowing iCalzetti et all (|2007[ ): SFR24 / (Mq yr^^) = 
1.27 X 10-38(L24/erg s'^-^^^. For NGC 1266, L24 = 
1.06 X lO^'^ erg s~^ on ce adjusted to our adopted dis- 
tance (|Temi et al.l 120091) . so the associated SFR24 « 1.5 
Mq yr~^ (the stellar continuum contribution to the 24^m 
emission is << 1%), within 33% of the SFRfir derived 
from the FIR emission. To measure the SFR from the 
PAH emission, we use the mid-IR fluxes from Falcon- 
Barroso et al. (in prep). We subtract off the stellar 
contribution from the 8.0/xm flux by scaling the 3.6//m 
image; the scale fa ctor of the 3.6/xm i mage can range 
from 0.22 to 0.29 (jCalzetti et "all 120071 ) with elliptical 
galaxies having typ ical values of 0.26 (jWu et al.l 120051 : 
IShapiro et al.l 12010 ). From the corrected PAH emission, 
we measure SFRpah ~ 0.57 Mq yr~^(Falc6n-Barroso et 
al. 2011, in prep) a factor of 4 lower than the SFRfir 
derived using the FIR emission. 

Comparing the SFR predicted from the molecular gas 
via the K-S relation to the SFRs derived from the 8^m, 
24/xm and FIR emission, we find that the predicted SFR 
is a factor of 1.4 to 5 times larger than the measured 
SFR, depending on which SFR tracer is used. However, 
NGC 1266 is well within the scatter of the K-S rela- 
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Figure 12. Chandra smoothed hard X-ray image truncated to 
only include 4-8 keV photons (greyscale and black contours), over- 
laid with the 5 GHz r adio continuum emission (cyan contours; 
IBaan fc Klockneil 120061 ) . The radio core and hard X-rays from 
Chandra are co-spatial within the expected positional accuracy of 
the Chandra data in the hard X-ray band (2"), convincingly iden- 
tifying NGC 1266 as harboring an AGN. 

tion (jLerov et al.l 120081 ). as well as the range of SFRs 
observed in the SAU RON and ATLAS^^ sample galaxies 
(jShapiro et al.|[20lol : Falcon-Barroso et al. in prep). 

The star formation-driven wind model of lMurrav et al.l 
(|2005h requires that the star formation rate be equal to 
or greater than the mass outflow rate for radiation-driven 
outflows, which we refer to as the Murray criterion. Our 
calculations for NGC 1266 show that even if we consider 
the higher SFR estimate (SFR « 2 Mq yr~^), our lower 
limit on the mass outflow rate (M « 13 Mq yr~^) still 
exceeds it by a factor of at least a few. This will only 
increase with more accurate accounting of the the total 
mass of the outflow. Therefore, the NGC 1266 molecular 
outflow does not fulfill the Murray criterion, and momen- 
tum coupling (i.e. radiation-driving) and star formation 
are unlikely to be the drivers of the outflow. 

4.3. AGN Feedback and Comparison to Other Molecular 

Outflows 

NGC 1266 is classified as a low-ionizati on nuclear 
emiss ion-line region (LINER) (_Moustakas fc Kennicuttl 
|2006( ). The 1.4 and 5 GHz radio continuum images reveal 
a prominent and unresolved point source in the nucleus 
(e.g. Figs [TT] and [T^ . C/iandra observations reveal hard 
X-ray emission centered on the radio point source, con- 
sistent with an obscured AGN (see Figure [T2|) . Given 
the overwhelming evidence for the presence of an AGN 
in NGC 1266, and that star formation is much less in- 
tense than in the few other galaxies exhibiting molecular 
outflows (see ij4.3p , an AGN-drivcn outflow seems a rea- 
sonable possibility. 

To investigate whether the AGN is any more 
capable, we use the relationship between the 



1.4 GHz l u minosi ty and total jet power from 
iBirzan et all (f2008l ): log(Pjet/10'*2 erg s"i) = 
0.351og(Pi.4GHz/102'' WHz-i) + 1.85. The total 
radio flux at 1.4 GHz is 9.3 x 10^" W Hz"\ rcsuhing 
in a total jet power of 6.1 x 10"*^ erg s~^. Comparing 
this to the Lmoch of the outflow of 1.3 x lO"'^ erg s~^, a 
coupling of 2% is required to drive the outflow with 
the AGN, a rather modest value. Therefore, the AGN 
does indeed generate sufficient energy to power the 
outflow through the mechanical work of the radio jet. 

There is a small sample of galaxies that have been 
shown to exhibit molecular outflows: M82, Arp 220, 
Mrk 231 and M51. Arp 220 is a starbursting major 
merger an d has a molecular gas surfac e density of 6 x 10^ 
Mq pc~^ (jMauersberger et 311119961 ). si milar to that of 
the rn olecular nucleus of NGC 1266. ISakamoto et al.l 
(|2009D detected an outflow that they estimate has a mass 
of 5 X 10^ Mq and a mass outflow rate of 100 Mq yr~^, 
with an outflow velocity of ~ 100 km s~^. These num- 
bers, though uncertain, are also comparable to what is 
observed in NGC 1266. Howe ver, the star forma tion rate 
in Arp 220 is > 200 Mq yr"! (jSoifer et al.lll987f ). almost 
two orders of magnitude greater than in NGC 1266, pro- 
viding a ready source of energy for the outflow and mean- 
ing that the Murray criterion is met. 

M82 has « 3 x 10^ Mq of molecular gas e ntrained in an 
outflo w with velocities up to « 230 km s^^ (jWalter et al.l 
120021 ). It is unclear whether this is sufficient for any gas 
to escape the galaxy. The mass outflow rate is « 30 
Mq yr~^, comparable to that estimated for the wind in 
NGC 1266, and rough ly equivalent to the global SFR of 
M82 ()Kennicuttlll998[ ). thus satisfying the Murray crite- 
rion. The outflows can thus be driven out via radiation 
p ressure. 

[Matsushita et all (j2007[ ) have shown recently that M51 
has entrained molecular gas in an outflow close to the 
AGN, with a total molecular mass of 6 x 10^ Mq, two 
orders of magnitude less than the mass estimate for the 
NGC 1266 outflow. The derived mass outflow rate of 
M51 is w 4 Mq yr~^. While this outflow rate is equiva- 
lent to the total SFR, the star formation in M51 is known 
to be distributed throughout the disk (jCalzetti et al.l 
120051 ). and thus is likely unable to effectively drive much 
of the nuclear mol ecular gas out. T he AGN luminosity of 
2 X lO'*^ erg s'^ (|Terashima fc Wilson. 2001 is also too 
low to sustain the molecular outflow via photon-driving. 
Therefore, M51 appears to be a scaled down version of 
the outflow wc observe in NGC 1266, requiring a similar 
mechanism, such as mechanical work done by a radio jet, 
to drive out the gas. 

Mrk 231, a nearby advanced major merger, shows 
a molecular outflow similar to that of NGC 1266 
(|Feruglio et aLllMol ). The Mrk 231 CO(l-O) spectrum 
exhibits a similar profile to that seen in NGC 1266, with 
broad wings requiring a nested Gaussian fit. This simi- 
larity seems to confirm that neither NGC 1266 nor Mrk 
231 is a unique case, and AGN feedback may be an ef- 
ficient means of removing th e molecular gas close to the 
AGN. iFeruglio etall ()2010f ) claim that Mrk 231 has a 
mass outflow rate of « 700 Mq yr~^ (assuming a density 
proflle of r~^), exceeding its star formation rate of 200 
Mq yr~-'^and thus requiring AGN feedback to be pow- 
ered. If the CO in the outflow is optically thick, as as- 
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sumed by iFeruglio et alj ()2010D , then the mass outflow 
rate in Mrk 231 is indeed too high for a star formation- 
driven wind. On the other hand, if the correct ij(CO)- 
to-M(H2) mass conversion for both the NGC 1266 and 
Mrk 231 outflows is the optically thin conversion (as as- 
sumed for NGC 1266), then A/(H2) for the Mrk 231 out- 
flow is a factor of w 5 lower, leading to a mass outflow 
rate of ^ 100 M© yr~^. This is sustainable using the 
Murray criterion, and is similar to Arp 220. Therefore 
the case for the molecular outflow in Mrk 231 requiring 
AGN feedback is uncertain. To obtain a robust estimate 
of the outflowing mass of Mrk 231, the opacity of the 
outflowing component must be determined. 

New stacked CO(l-O) observations of local ultra- 
lumin ous infrared galaxies (ULIRGs) by iChung et al.l 
(|2011[) also show clear evidence that high- velocity broad 
emission at up to 1000 km s~^ from the systemic velocity 
is a common feature in these systems. In particular, the 
authors find that only the starburst-dominated ULIRGs 
exhibit broad wings, while AGN~dominated ULIRGS 
(identified by their optical spectra) do not. These results 
are thus clear evidence that in most observed molecular 
outflows, photon-coupling with the star formation is the 
principle driver. 

While the aforementioned galaxies all exhibit molec- 
ular outflows, M51 is the only other bonafide case of a 
non-photon-driven molecular outflow. However, it has 
almost two orders of magnitude less molecular gas being 
driven out than NGC 1266. Mrk 231 might well be an 
equivalent system, but the mechanisms capable of driv- 
ing the gas out depends on the CO luminosity-to-H2 
mass conversion in the outflow. The most acute con- 
trast between NGC 1266 and the other molecular out- 
flow sources is its lack of interaction. All other galaxies 
are undergoing interactions; major m ergers in the case o f 
Arp 220, Mrk 231 and some of the (jChung et al.ll20lH) 
ULIRGs. NGC 1266 shows no evidence of having un- 
dergone an interaction at any time in the recent past, 
with unperturbed stellar kinematics, no companion, and 
a lack of H I emission internal or external to the galaxy. 
This means that while we can point to an interaction as 
the triggering mechanism in the other known molecular 
outflows, the striking absence of such a perturbing event 
makes NGC 1266 even more remarkable. 

5. SUMMARY AND CONCLUSIONS 

We report the detection of a massive centrally- 
concentrated molecular component (Mcvc = 1.1 x 10^ 
Mq) and a powerful molecular outflow (Mh2, outflow = 
2.4 X 10^ Mq) in the field SO galaxy NGC 1266. The 
maximum velocity of the wind exceeds the escape veloc- 
ity so that at least some gas will escape the galaxy to 
energize the IGM, and the total neutral gas mass out- 
flow rate AI w 13 Mq yr~^. The star formation rate in 
NGC 1266 appears to be incapable of driving this out- 
flow; the AGN appears to be the main driving mecha- 
nism, as for H I outflows observed in radio galaxies. It is 
however unclear how the gas lost its angular momentum 
to fall so completely into the nucleus. 

The central molecular gas contains a rapidly rotating 
nuclear disk of « 60 pc radius enshrouded in a diffuse 
molecular envelope, and the molecular outflow emerges 
normal to the disk plane. Although the molecular nu- 



cleus is compact and very near the AGN, multiple es- 
timates of the star formation rate point to a low star 
formation efficiency, but the uncertainties in these SFRs 
place NGC 1266 within the scatter of the K-S relation 
(assuming that LpiR is not dominated by the AGN). 

If the gas in the nucleus is the source of the molecu- 
lar outfiow, the gas depletion timescale is 85 Myr, short 
enough for all ETGs to go through such a phase while 
remaining consistent with having found only one such 
case among the 260 ETGs in the ATLAS^^ sample. 
NGC 1266 is the brightest detection, however, so there 
might be more outfiows in the sample, undetected so far 
because of lower S/N ratios. 

While molecular outflows have been identified in a 
handful of galaxies, NGC 1266 is the only one that shows 
no evidence of having undergone an interaction, leaving 
the mechanism to transport the molecular gas into the 
center unknown. 
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